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Abstract: We study the determination of super symmetric parameters at the LHC from a 
global fit including cross sections and edges of kinematic distributions. For illustration, we 
focus on a minimal supergravity scenario and discuss how well it can be constrained at the 
LHC operating at 7 and 14 TeV collision energy, respectively. We find that the inclusion 
of cross sections greatly improves the accuracy of the SUSY parameter determination, and 
allows to reliably extract model parameters even in the initial phase of LHC data taking 
with 7 TeV collision energy and 1 fb _1 integrated luminosity. Moreover, cross section 
information may be essential to study more general scenarios, such as those with non- 
universal gaugino masses, and distinguish them from minimal, universal, models. 
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1. Introduction 

Experiments at the Large Hadron Collider (LHC) will, for the first time, directly explore 
supersymmetry (SUSY) [1, 2] at the Terascale. Provided that evidence for supersymmetry 
has been established, a major challenge will be to determine the Lagrangian parameters of 
the theory, such as the SUSY particle masses, their spins and couplings. These TeV-scale 
parameters provide essential information on the scheme of supersymmetry breaking and 
need to be determined with the highest possible accuracy. 

A generic LHC signature for super symmetric models with R-parity conservation [3, 4] 
is that of cascade decays of heavy squarks and gluinos which terminate in the lightest 
super symmetric particle (LSP). In many SUSY models the LSP is weakly interacting and 
provides a viable dark matter candidate [5, 6]. However, a weakly interacting LSP escapes 
detection and thus results in missing energy signatures. It is a considerable challenge to 
reconstruct the sparticle momenta in such cascade decays with missing energy at the LHC 
and to determine sparticle masses and quantum numbers. The standard technique for 
analyzing SUSY cascade decays is to consider invariant mass distributions of the final- 
state leptons and jets, see e.g. [7, 8, 9, 10]. The kinematic endpoints of these distributions 
are fixed by the masses of the sparticles in the decay chain and yield model-independent 
information on part of the SUSY mass spectrum. The endpoints can be used as input for 
global SUSY parameter fits of LHC data [11, 12, 13, 14, 15] which determine the high- 
scale model parameters and thereby test mechanisms of SUSY breaking. (See [16] for a 
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discussion of various fitting methods and tools). However, the endpoints of distributions 
are only sensitive to the kinematics and can not constrain important SUSY parameters 
such as tan/3 very well. Moreover, as the endpoints are determined by the difference of 
masses, there are ambiguities, i.e. points in the SUSY parameter space that have different 
spectra but similar endpoints. There are further ambiguities from cascade decays with 
identical final states, but different intermediate particles. It is thus of vital importance to 
consider additional information for the determination of SUSY parameters at the LHC. 

In this paper we address the impact of event rates, i.e. cross sections and branching 
ratios, on SUSY parameter fits of LHC data. The cross sections are highly sensitive to 
the masses of the squarks and gluinos that trigger the cascade decays. Branching ratios, 
on the other hand, depend sensitively on the masses and mixings of SUSY particles fur- 
ther down the decay chain. The inclusion of event rates is technically challenging as the 
multi-dimensional parameter scan in SUSY fits requires computationally extremely fast 
and reliable theoretical estimates for cross sections, branching ratios and the effect of ex- 
perimental cuts. The most straightforward and most flexible approach of estimating event 
rates from a Monte Carlo simulation for every point in the parameter fit is prohibitively 
slow [11], so that new strategies need to be developed. Moreover, Monte Carlo calculations 
are generally based on leading-order perturbation theory, resulting in a substantial theo- 
retical uncertainty in the prediction of event rates. Finally, the statistical fluctuations of 
Monte Carlo estimates can cause oscillations during the ^-minimization in gradient based 
global fits and lead to unstable results [15]. 

We propose a new approach to include event rates into SUSY parameter fits. Our 
method is based on a simple parametrization of cross sections and acceptances 1 , and does 
not involve Monte Carlo simulations during the ^-minimization. It is thus fast and re- 
producible, and it incorporates state-of-the-art higher-order cross section predictions with 
small theoretical uncertainties. Branching ratios, on the other hand, can be evaluated very 
quickly with existing computational tools (see e.g. [17]). The information on event rates, 
i.e. cross sections x branching ratios x cut acceptances, has been implemented as an addi- 
tional observable into Fittino [18] and will become part of the next official release of the 
Fittino program package [19]. 

We expect event rate information to be particularly valuable for general SUSY sce- 
narios beyond minimal supergravity, for example those which involve three-body decay 
modes, so that mass reconstruction via the standard kinematic endpoints is difficult, see 
e.g. [20, 21]. Moreover, parameter determinations of theories with a general, non-universal, 
gaugino mass pattern may play a crucial role to distinguish models of supersymmetry break- 
ing [22]. As we shall see, cross sections and branching ratios do not only add information 
to stabilize a fit with the larger set of parameters in non- universal models, they also allow 
to determine specific gaugino masses such as M3 much more reliably. Finally, production 
cross sections are essential to distinguish supersymmetric theories from alternative new 
physics models, such as universal extra dimensions [23], which predict similar cascade de- 

x We use the terms "acceptance" or "cut acceptance" to denote the fraction of events that passes the 
experimental selection cuts. Our acceptance estimates are based on parton level calculations and do not 
yet include detector effects. 
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cay signatures at the LHC, but involve new particles with spin quantum numbers different 
from the SUSY particle spectrum [24]. A universal extra dimension model could not be 
distinguished from a supersymmetric theory based on the kinematic information from end- 
points of distributions, but would have dramatically different cross sections because of the 
different particle spins [25, 26, 27]. 

To exemplify the impact of event rates on SUSY parameter fits we have considered the 
standard SPSla minimal supergravity (mSUGRA) scenario [28], which has been studied 
in great detail in previous analyses. SPSla provides a rich signature at the LHC and can 
thus be constrained well by standard measurements of kinematic endpoints. It is thus 
not a scenario where cross sections and branching ratios are expected to have the largest 
impact. On the other hand, SPSla allows for a quantitative analysis of the effects of 
event rates on the parameter fit with realistic error estimates, as discussed in detail in 
section 3. Due to the relatively light spectrum it is also one of the SUSY benchmark 
points that can be studied in the initial phase of the LHC with 7 TeV collision energy and 
low integrated luminosity. Specifically, we investigate the standard signatures of SPS la- 
type models including two or more jets, missing transverse momentum and two leptons 
of the same flavour but opposite sign. Details of the cuts that specify the signatures 
will be given in section 2. We use Fittino to determine the parameters of the SPSla 
mSUGRA scenario from the measurement of kinematic edges and event rates at the LHC, 
considering 7 TeV collision energy with 1 fb _1 integrated luminosity and 14 TeV with 
both 1 fb _1 and 10 fb _1 . As demonstrated in section 3 the inclusion of rates in general 
stabilizes the fit and significantly improves the error on the mSUGRA parameters, in 
particular the universal gaugino mass Mi/% and the ratio of vacuum expectation values 
tan f3. Moreover, we find that the inclusion of rates is crucial for the determination of 
mSUGRA parameters in the initial phase of LHC data taking with 7 TeV collision energy 
and 1 fb _1 integrated luminosity. We also consider a fit where we determine the gaugino 
mass parameters Mi, M2, M3 individually at the GUT scale, instead of M 1 / 2 , and show 
that rate information is important to improve the accuracy of the parameter determination 
for more general, non-universal, models. 

The paper is structured as follows. Section 2 presents details of our method to esti- 
mate and parametrize the cross sections and the effects of experimental cuts. The SUSY 
parameter fits are discussed in section 3. We conclude in section 4. 

2. Cross sections and cut acceptances 

The prediction of rates involves the calculation of the production cross section for squarks 
and gluinos, which dominate the inclusive SUSY signal at hadron colliders, the branching 
ratios of the supersymmetric particles in the decay chains, and an estimate of the effect of 
a certain set of experimental cuts. We have considered two powerful and widely studied 
SUSY signatures [7, 8, 9, 10]: 

- the inclusive signal of two or more hard jets with each PTjet > 50 GeV, \f]\ < 2.5 as 
well as missing transverse energy of Ifi T > 100 GeV; 
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- the exclusive signal of two opposite-sign same- flavour leptons (e or fi), each satisfying 
PT,e > 10 GeV and \rj\ < 2.5, combined with the above signal of two or more hard jets 
and missing transverse energy. We remove the background from the leptonic decays 
of tau leptons, charginos and bosons in the standard way by subtracting events 
with opposite-sign different-flavour lepton pairs, see [7, 8, 9, 10]. 

Assuming a typical mSUGRA-like mass spectrum with m g > rriq, the hard jets result from 
two-body decays of gluinos and squarks, g — >■ qq and q — > qx, respectively, where x denotes 
a chargino or neutralino, i.e. a model-dependent linear combination of the charged and 
neutral gauginos and higgsinos. The leptons are produced in chargino/neutralino decays 
further down the cascade chain, e.g. x% ~^ $R ^Xi> as illustrated in Fig. 1. 

In the SPSla scenario which we con- 
sider for illustration in this paper, the 

neutralinos X3 4 are mainly higgsino, so >- - 

that the decays q — > 17X3,4 are strongly 1 
suppressed. Decays into on-shell Z or 
Higgs, x° -> z Xi or X2 -> h Xi, need a 
sufficient x\ ~ X? mass splitting and are 
thus not accessible for many mSUGRA Fi g Ure \. The cascade decay of a squark leading 
scenarios, including SPSla. Such decays to the leptonic signal, j denotes the quark which 
into Z or Higgs will be considered in the should lead to a hard jet, n denotes the near lepton, 
future, together with other decay modes / denotes the far lepton. 
with potential significance for more gen- 
eral SUSY models, but are not relevant for the results shown in this paper. 

The calculation of the event rates proceeds in various stages. First, for each point in 
the SUSY parameter space, the mass spectrum and branching ratios are calculated by a 
spectrum generator. The masses and branching ratios are highly model dependent, but 
they can be evaluated quickly and can thus be calculated during the ^-minimization for 
every point in the parameter space. Our results are based on SPheno [29], which is the 
default generator used by Fittino. Given the mass hierarchy, we first check which cascade 
decays can contribute to the signal. The corresponding event rates are proportional to 
the production cross sections for squarks and gluinos. These cross sections essentially only 
depend on the squark and gluino masses and on no other SUSY parameters or model 
assumptions. They can thus be calculated once and for all, including all available higher- 
order QCD corrections, stored in a (nig, rrig)-gx\d and read out quickly during the fit, see 
section 2.1. We calculate the cut acceptances using a combination of numerical results 
stored in look-up tables and analytical calculations of jet energies and lepton energies in 
the squark rest frame. The calculation of the cut acceptances is less straightforward, and 
details will be presented in section 2.2. For each contributing cascade, the cross-section 
for the production of the relevant colored sparticles is multiplied by the relevant branching 
ratios and the cut acceptance, and the total event rate is passed on to Fittino as an 
observable entering the SUSY parameter fit. 

In the following, we will describe the calculation of cross sections and cut acceptances 
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in some detail. 

2.1 Cross sections 

In supersymmetric models with R-parity conservation, squarks and gluinos are produced 
in pairs or associated pairs at hadron colliders, pp — > qq, qq*,qg, gg + X. Here we suppress 
the SU(2) quantum numbers of the squarks q = (qx, q~R) and do not explicitly state the 
charge-conjugated processes. The squark and gluino production cross sections are known 
including next-to-leading order (NLO) SUSY-QCD corrections [30, 31, 32, 33], the sum- 
mation of soft gluon emission [34, 35, 36, 37, 38], as well as electroweak contributions and 
corrections [39, 40, 41, 42, 43, 44, 45]. The strength of the SUSY-QCD interactions and 
thus the production rate is set by the gauge and the Yukawa couplings of the qqg and qqg 
interactions, respectively. The two couplings are required by supersymmetry to be equal, 
so that the LO and NLO QCD squark and gluino parton cross sections can be predicted 
unambiguously in terms of the QCD coupling g s and the squark and gluino masses, with- 
out any further SUSY model dependence. The NLO-SUSY QCD corrections have been 
implemented in the public computer code Prospino [46] and form the basis of our cross 
section prediction. Note that the strong Yukawa coupling between top quarks, top squarks 
and Higgs fields gives rise to potentially large mixing effects and mass splitting in the top 
squark sector. The mixing angle enters the top-squark cross section at NLO, however, the 
dependence is numerically very weak [33]. Soft gluon resummation leads to a reduction 
of the scale dependence of the cross section prediction and an enhancement of the NLO 
QCD cross section for heavy squarks and gluinos with masses fh > 1 TeV [38]. These 
effects will be included in future studies, but are not essential for the results shown in 
this paper. Electroweak contributions and corrections introduce a dependence on further 
supersymmetric parameters. Their impact on the inclusive cross sections summed over all 
squark species is generically small. However, interference effects between the exchange of 
electroweak gauge bosons and QCD contributions can become large for the production of 
two SU(2) doublet squarks pp — > (Ilq'l [40, 41]. For SPSla-type scenarios, the electroweak 
effects are moderate, but they may be larger for more general models, in particular those 
without gaugino mass unification. Given the residual QCD uncertainties of the produc- 
tion cross section discussed below, we can neglect electroweak contributions for the study 
of SPSla presented in this paper. However, such electroweak effects will be addressed in 
future extensions of this work. 

The calculation of the squark and gluino production cross sections involves the cal- 
culation of the corresponding parton cross sections, including higher-order SUSY-QCD 
corrections, and the convolution with parton distribution functions. Numerically, this is a 
time-consuming task, even at leading order, and would make the inclusion of cross sections 
in SUSY parameter fits prohibitively slow. Hence, we have calculated the cross sections as 
a function of rriq and m g and stored the results in a grid. As the cross sections are smooth 
functions of the sparticle masses, the values between grid points can be interpolated reli- 
ably. When considering the leptonic signature, we have to distinguish the SU(2) quantum 
numbers of the squark, as decays from q~R and q~i are in general very different. We have thus 
calculated the LO cross sections for the production of squarks with definite SU(2) quantum 
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number, qL and q~R, averaging over the masses of q G (u, d, c, s) for each SU(2) quantum 
number. Top and, to a lesser extend, bottom squarks mix to form mass eigenstates and 
are thus treated separately, as mentioned above and explained in detail in Ref. [33]. The 
NLO SUSY-QCD corrections are taken into account through If -factors, K = ctnlo/^lO) 
as provided by Prospino. Note that the calculations [30, 31, 32] implemented in Prospino 
sum over squark SU(2) quantum numbers and do not provide separate if-factors for the 
production of <jl and qR. We thus assume that the if-factors do not depend significantly on 
the squark SU(2) quantum number. Furthermore, we average the if-factors for the qq and 
qq* channels. For our numerical results, we have set the renormalization and factorization 
scales to the average mass of the produced sparticles and adopted the 2008 MSTW parton 
distribution functions [47]. 

To illustrate the sensitivity of the SUSY cross section to the squark and gluino masses, 
we show in Fig. 2 the NLO cross sections for pp — > qq,qq* ,qg,gg + X at the LHC with 
14 TeV collision energy. In the Figure we sum over q G (u,d,c,s,b) and L/R quantum 
numbers, and average over the squark masses. 
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Figure 2: NLO QCD cross-section for inclusive squark and gluino production at the LHC (14 TeV) 
in pb, as a function of the gluino and average squark masses. 

Varying (m^, rrig) in the range between 200 GeV and 2 TeV, the cross section changes 
by seven orders of magnitude. It is evident that the sensitivity of the cross section to the 
sparticle masses should play an important role in SUSY parameter fits. 

Let us finally comment on the theoretical error of the cross section prediction. The 
renormalization and factorization scale uncertainty of the NLO QCD cross section is < 10% 



- 6 - 



for squark and gluino masses below approximately 1 TeV [32]. This uncertainty could be 
reduced further by taking into account soft gluon resummation [38]. In addition, there is 
the uncertainty due to the parton distribution functions, which, however, is estimated to 
be below 5% for sparticle masses less than about 1 TeV, see e.g. [47, 48, 49]. We thus 
assume an overall theoretical uncertainty of 15% on our NLO cross section prediction. We 
shall discuss in section 3.1 how this uncertainty enters the SUSY parameter fit. 

2.2 Cut acceptances 

The possibility of parameterizing the cut acceptances in a generic way hinges upon fac- 
torizing the production of the colored sparticles and the subsequent cascade decay. This 
factorization relies, of course, on the narrow width approximation, which is appropriate 
as the width of the SUSY particles is in general small, T <C m? We can thus break the 
problem of estimating cuts for the full production % decay process down to calculating the 
decay distributions for jets and leptons analytically in the squark rest frame and to nu- 
merically estimating the effect of boosting the particles from the squark rest frame to the 
laboratory frame. The distributions of the squark decay products depend on the SUSY 
scenario and on many SUSY parameters, but they can be calculated analytically, and thus 
evaluated quickly during the SUSY parameter scan. The impact of the boost from the 
squark rest frame to the laboratory frame depends on the dynamics of the production 
process and is in general difficult to obtain analytically. However, this part can be treated 
numerically and stored in a look-up table, as it only depends on the production dynamics 
and thus the squark and gluino masses, and not on any specific details of the SUSY sce- 
nario. Our method will be outlined briefly in sections 2.2.1 and 2.2.2; more details can be 
found in [51]. Currently, the parametrization of cut acceptances is implemented for sig- 
nals that consist of sequences of two-body decays with the usual mSUGRA mass hierarchy 
rrig > rriq > m^a > > m^o. More general scenarios will be considered in future work. 

Even though the inclusive cross-sections for the production of colored sparticles have 
large NLO QCD corrections in general, the NLO effects on the differential distributions and 
thus on the typical boosts of the colored sparticles are small [32]. Hence, for the estimate 
of the acceptance we rely on leading order calculations. 

2.2.1 Jet and missing energy cuts 

Unfortunately the missing energy cut is very difficult to calculate analytically for a generic 
decay cascade, and we resorted to a numerical grid of acceptances. The grid has been 
obtained by a simple parton-level Monte Carlo simulation, decaying all particles by phase 
space and ignoring spin correlations. This is a legitimate approximation, as we average over 
charges in the final state. We find that the effects of intermediate decays from the squark 
to the lightest neutralino tend to average out, and the missing energy cut acceptance is 
approximated well as a function of the hard process (giving the typical squark boosts) and 
just the mass difference between the squark and the lightest neutralino. The acceptance 
grid can thus be parametrized by three masses only: rrig, rrig and rn^o. 

2 See Ref. [50] for new physics scenarios where finite width effects are important; these scenarios are, 
however, not relevant for the type of models studied in this paper. 
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Since all the jets from the decays of the gluinos and squarks were simulated in the 
process of calculating the missing energy cut acceptances, the jet cut acceptances were also 
taken from these simulations, though in principle the jet cuts could have been dealt with 
in the same manner as the lepton cuts described below. It is intended to implement this 
in the future, to allow for the jet cuts to be specified by the user, rather than hard-coded. 

2.2.2 Lepton cuts 

As argued above, we can estimate the impact of lepton cuts on the production <g) decay 
process by calculating the decay distributions to leptons analytically in the squark rest 
frame and by numerically estimating the effect of the squark boost to the laboratory frame. 
We have derived analytical expressions for the distributions of near and far leptons in the 
squark rest frame, which can be evaluated very quickly during the global fit for every 
point in parameter space. To estimate the lepton acceptance in the laboratory frame, 
we numerically calculate the acceptance for a generic massless lepton at a given energy 
in the squark rest frame as a function of the lepton energy and perform a boost to the 
laboratory frame by means of a simple Monte Carlo calculation. Since the squark is a spin- 
zero particle, the momentum of any one of its final decay products has a flat solid angular 
distribution in the squark rest frame, though of course there are correlations between the 
momenta of the decay products. However, we find that the correlations are weak, and that 
merely multiplying together acceptances for each jet or lepton produces reasonable overall 
acceptances. The numerical estimate of boosting to the laboratory frame depends on the 
dynamics of the squark and gluino production process only and is stored in a grid as a 
function of the squark and gluino masses and the lepton energy. We finally multiply the 
generic acceptance estimate with the number of leptons at a given energy, as obtained from 
the analytic calculation of the squark-to- lepton decay. Though the generic acceptances were 
obtained for a particular choice of transverse momentum cut (10 GeV), they scale simply 
and can be re-used for different choices of the cut. 

2.2.3 Verification 

To validate our acceptance cut estimates, we performed full parton-level simulations, inclu- 
ding spin correlations, with Herwig++ [52] for a set of mSUGRA points chosen randomly 
with flat priors, restricted then to those points with spectra such that the gluino decayed to 
on-shell squarks, and that the lightest neutralino was the LSP. The jet and missing energy 
acceptances from the Herwig++ simulations were compared to those calculated using the 
methods described above, and were found to agree within 5% or better. We thus attach 
an uncertainty to our acceptance cut estimate of 5%. 

3. Numerical results 

The calculation of the event rates described in section 2 is available in the form of a C++ 
code and has been incorporated into a currently private modification of Fittino. It will 
become part of the next official release of the Fittino program library. Fittino is a 
program which attempts to find the best fit for Lagrangian parameters of super symmetric 
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models given a set of observables with uncertainties and their covariance matrix. One can 
choose to perform a fit using simulated annealing or Markov chains (see [15] and references 
therein), and one can choose to fit to various sets of Lagrangian parameters defined at 
either the TeV scale or the GUT scale. Since the results presented below were obtained 
by performing fits with Markov chains to parameters defined at the GUT scale, we briefly 
outline this process. 

Firstly, Fittino reads in its input file, in which all the observables to be used are 
defined, with their nominal values, uncertainties and correlations. These observables, and 
the corresponding values adopted in our analysis, are presented in section 3.1. Fittino 
also reads in the parameters which should be fitted, and a set of starting values for these 
parameters. Since our aim is to show how the addition of event rates as observables helps 
to reduce uncertainties on the parameters of the fit, we began the Markov chains at the 
point which should be found, rather than some other point in the parameter space. At 
this stage, our code is initialized, and loads the cross-section and acceptance grids into 
memory. Secondly, Fittino begins the process of scanning the parameter space. At each 
step, it selects a random point near its current point. Then SPheno [29] is used to calculate 
the TeV-scale spectrum and branching ratios. Fittino uses this spectrum to calculate 
theory predictions for the observables (either with internal code, including our extension 
to rates, or by calling external programs). These predictions are then compared to the 
input values and uncertainties to calculate a likelihood for this new point, which is then 
used by a Metropolis algorithm [53] to decide whether Fittino moves to this point or 
rejects it and stays where it is. Finally, the Markov chain is analyzed to calculate a best 
fit for the parameters and their uncertainties. 

3.1 Fit inputs 

Fittino allows for the inclusion of a wide variety of low-energy and hadron as well as lep- 
ton collider observables. For the case study presented here, we restrict ourselves to those 
observables which are expected to be measured at the LHC. We consider the minimal su- 
pergravity SPSla benchmark point, with parameters Mq = 100 GeV, M 1 / 2 = 250 GeV, 
Aq = —100 GeV, tan/3 = 10 and sign(/x) = +1 [28], corresponding, e.g., to the TeV- 
scale masses m~ g = 606 GeV, rriu L = 559 GeV, m^ L = 177 GeV, m^o = 181 GeV and 
m^o = 97.1 GeV. SPSla has been studied in great detail in experimental simulations by 
the ATLAS collaboration [10]. Because of the comparably light spectrum, this benchmark 
point provides a rich phenomenology, even at low LHC collision energy and luminosity. 
Furthermore, SPSla has been adopted as an input for previous SUSY parameter determi- 
nations from LHC data [13, 15], the results of which may be compared to our's. 

The standard set of LHC observables for SUSY mass determination comprises the end- 
points of the four invariant mass distributions that can be constructed from the jet and 
the two leptons in the cascade depicted in Fig. 1: 

- m^ ax , the dilepton invariant mass edge, 

- m™$ x , the jet-dilepton invariant mass edge, 

- mj^ w , the jet-lepton low invariant mass edge, and 
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- rn^f 1 , the jet-lepton high invariant mass edge. 

The definition of these edges and their relation to the sparticle masses in the decay chain can 
be found in [8]. The m^ ax , m^ 1 , and rn^ h endpoints have been analyzed in detail by 
both the ATLAS and CMS collaborations for different SUSY scenarios [10, 54, 55, 56, 57]. 
At the SPSla benchmark point, they can be measured at the LHC with a high accuracy 
of better than 5%, even at a low luminosity of 1 fb _1 [10, 15]. The estimated statistical 
uncertainties on the measurements of the edges are collected in Table 1 for different LHC 
collision energies and luminosities. We will present results for SUSY parameter fits that 
involve the four standard edges (labeled group I) and the rates for LHC data analysis at 
7 TeV collision energy and 1 fb _1 integrated luminosity, and at 14 TeV for both 1 fb _1 
and 10 fb -1 . 

For 14 TeV collision energy, we also consider a set of observables which involve third- 
generation particles, the lower endpoint of m q u and the stransverse mass [58, 59] of the 
right-handed up squark: 

- m^', the jet-dilepton threshold invariant mass edge, 

- m^, the squark stransverse mass, 

- m™ T ax , the di-tau invariant mass edge, 

- m™ b , the weighted top-bottom invariant mass edge, and 

- rjj fBR , the ratio of selectron- to stau-mediated decays. 

This second group of observables (group II) is defined in [10, 15] and is much more challeng- 
ing experimentally than the edges of group I. We are only aware of an ATLAS study [10, 15] 
which quantifies the experimental accuracy to be expected at the LHC operating at 14 TeV. 
In any case, also this second group of measurements may contribute to the SUSY parameter 
determination at the highest LHC energy and will be included in our analysis below. 

The third group (group III), finally, is a set of observables which may only be measur- 
able with at least 10 fb -1 of data at 14 TeV, and comprises [10, 15] 

- Am^o, the mass difference between the gluino and the LSP, 

- m™^ u , the dilepton invariant mass edge from the decay of & xl, 

- m$g-, the 6-tagged jet-dilepton threshold invariant mass edge, and 

- m/j, the mass of the lightest neutral scalar Higgs boson. 

For all SUSY parameter fits we include the two types of event rates defined in section 2 
as additional observables: 

- Rjj$ T , the inclusive event rate for at least two hard jets with missing transverse en- 
ergy, and 

- R££jjp T , the exclusive event rate for at least two hard jets with missing transverse 
energy plus a pair of opposite-sign same-flavour light leptons. 
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observable 



nominal statistical uncertainty 

value for 7 TeV/1 ftr 1 for 14 TeV/1 ftr 1 for 14 TeV/10 ftr 1 



group I 

ruff* 80.4 4.4 1.5 0.43 

452.1 36.0 12.0 3.6 

mjf 318.6 19.7 6.5 3.0 

mjf h 396.0 13.5 4.5 3.9 

group II 

mjg- 215.6 - 22.8 4.1 

m| 2 531.0 - 16.9 5.3 

m™ ax 83.4 - 10.8 3.4 

mf b 359.5 - 37.0 11.7 

r ifBR 0.076 - 0.008 0.003 

group III 

Am s o 507.7 - - 11.8 

!/ai 

m™^„„ 280.6 - - 10.8 

m$$- 195.9 - - 17.0 

m h 109.6 - - 1.2 



Event rate [fb] 7 TeV 14 TeV 

nominal value uncertainty nominal value uncertainty 
R mT 4.6 xlO 3 9.1 xlO 2 4.8 xlO 4 9.5 xlO 3 

Rttjjpr 1-6 xlO 2 3.2 xlO 1 1.5 xlO 3 3.0 xlO 2 



Table 1: LHC observables for SPSla used as inputs to the SUSY parameter fits. The nominal 
values for masses and branching ratios have been obtained with SPheno. The uncertainty estimates 
on the observables of groups I, II and III are based on [15] and have been rescaled as described in the 
main text. Note that the uncertainty estimates for 7 TeV are not based on a detailed experimental 
simulation but on a simple extrapolation from experimental studies at 14 TeV. The event rates 
include the NLO squark and gluino production cross sections, the branching ratios and the cut 
acceptances, as described in sections 2.1 and 2.2, respectively. The overall uncertainty on the rates 
is assumed to be 20%, see section 2.1 and the main text below. Dimensionful quantities are given 
in units of GeV for masses and invariant mass endpoints, or fb for event rates. 

The theoretical uncertainty of our rate prediction includes 15% uncertainty on the NLO 
calculation of the squark and gluino cross section and 5% uncertainty on our acceptance es- 
timate. Further theoretical uncertainties arise from the renormalization group running [60], 
the conversion of Lagrangian parameters into physical masses [61] and the calculation of 
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the branching ratios, all obtained at finite order in perturbation theory. These effects are, 
however, estimated to be small compared to the 15% production cross section uncertainty 
and can safely be neglected. For SPSla, the statistical uncertainty on the rate is also sig- 
nificantly smaller than the theory uncertainty: even at 7 TeV collision energy with 1 fb _1 
integrated luminosity, we expect about 150 signal events for the exclusive signature includ- 
ing two opposite-sign same-flavour leptons. We thus simply assume a conservative overall 
rate uncertainty of 20% for our analysis of SPSla. 

All observables are collected in Table 1 together with their estimated statistical uncer- 
tainties. The overall experimental uncertainties that enter the fit also include a systematic 
uncertainty given by the jet and lepton energy scale, which is assumed to be 5% (1%) 
and 0.2% (0.1%), respectively, for 1 (10) fb _1 [10, 15]. The uncertainties on the endpoints 
related to the lepton and jet energy scale are considered 100% correlated among different 
measurements. Further systematic uncertainties for some of the observables of groups II 
and III can be found in [15]. The statistical uncertainties in Table 1 are estimated by 
rescaling the uncertainties for 14 TeV collision energy listed in [15] by \/Rn, where Rn is 
the ratio of expected numbers of events passing appropriate cuts when going from LO at 
14 TeV (on which the simulations in [15] are based) to NLO at 14 TeV or NLO at 7 TeV. 

Our fits are obtained from Markov chains with 10 5 steps. The frequentist interpretation 
of these Markov chains is used to calculate the errors of the parameters, see [15] for details. 
Including rates slows down the computation by only about 30%; further optimization is 
possible and will be implemented in the public release of the code. 

3.2 Universal mSUGRA fits 

In this section and section 3.3 below we present the main results of this paper and demon- 
strate that rates significantly improve the fits of SUSY parameters at the LHC. We first 
discuss the determination of the universal mSUGRA parameters Mq, My 2 , tan/3 and Aq 
(with sign(/x) = +1 fixed) for the SPSla benchmark scenario. We show results for 7 TeV 
LHC collision energy with 1 fb _1 luminosity, and 14 TeV with both 1 fb _1 and 10 fb -1 . 

Let us first consider LHC data taken at 7 TeV and 1 fb -1 . As mentioned above, 
it is not clear which observables beyond the four standard edges can be measured with 
what accuracy during such an initial phase of LHC data analysis. We thus address the 
question of what information can be obtained from the four edges and the event rates only. 
Remarkably, we find that we can constrain the universal scalar and gaugino masses with an 
error of 10% and 3%, respectively. The results of our fit are collected in Table 2. Moreover, 
the fit reproduces the correct value of tan/3 with a reasonable error. Note that the tri- 
linear coupling Aq remains essentially unconstrained. The impact of the rate information is 
crucial: including the four edges of group I alone leads to very large errors and an unreliable 
fit, as one can judge from the profile likelihood contours shown in Figure 3. (Note that 
we needed 2 million Markov chain steps in total to derive the plot showing the fit without 
rates, while the fit with rates has been obtained with 10 5 steps only.) 

At 14 TeV collision energy and 1 fb" 1 luminosity, the statistical uncertainty on the 
measurement of edges is considerably reduced, see Table 1. We thus obtain a better 
determination of the mSUGRA parameters with rates and the basic edges of group I, in 
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Mq [GeVJ 


Mi/2 [OeVJ 


tan/3 


A [CjreVJ 


SPSla 


100 


250 


10 


-100 


7 TeV and 1 fb -1 










1 1 T*P"4"^G 

x latco 


qq n + 9 - 9 
yy.u _g j 


ZOU.U _ 6 5 


10 7 + 4 - 
iu -' -8.8 


crcr 9 +1048 
00 - z -254 


14 TeV and 1 fb" 1 










I + rates 
I + II, xatetf 
I + II + rates 


99.7 

99.8 tH 
99.8 ±££ 


251.1 
249.7 tH 
251.3 tH 


n.2 tn 

io.i ±1:1 

10.7 ±3-1 


— 50 Q + 1233 
ou - y -350 

_q 4 1 +1610 
y4 -! -216 

00 -' -233 


14 TeV and 10 fb" 1 










I + rates 
I + II, xates 
I + II + rates 


100.0 ±|;1 

100.1 t\- 7 9 

100.3 ±J;jj 


250.7 ±|;g 
250.4 ±{;f 
250.4 ±*;| 


ii.o tli 

10 1 + 11 

-LU- 1 -1.0 

io.2 tli 


-63.3 till 
"89.8 ±™ 

-96.5 tin 


I + II + III, -rates' 
I + II + III + rates 


100.2 tH 
100.1 ±J;j? 


250.3 ±J;£ 
250.3 ±J;£ 


io.i ±°:i 

10.3 t% 


"94.6 ±1^ 

-90.3 tir 7 



Table 2: Fits to universal mSUGRA parameters for SPSla, with ("+rates") and without ("xateS") 
event rates as an observable. The symbols I, II and III refer to the inclusion of the groups of 
previously considered observables (mainly edges) defined in the main text. 




Figure 3: A% 2 = — 21n£ + 21n£ max contours showing M against -Mi/ 2 for 7 TeV/1 fb _1 data. 
Fits are based on the four standard edges of group I without rates ( "I, xateS" , left) and with rates 
("I + rates", right). C is the two-dimensional profile likelihood and £ max the global maximum of 
the likelihood. The black dotted contours represent A% 2 = 1 contours. Sec [15] for more details. 



particular for the common scalar mass Mq. Note that the errors on tan (3 and M\ / 2 in the fit 
"I + rates" in Table 2 would be twice as large had we based our rate estimates on LO cross 
sections with the corresponding 100% uncertainty. Adding the observables of group II, 
which involve information on third-generation particles, the lower endpoint of m q u and 
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I + rates 




M„ 2 (GeV) 



I + II, xates" I + II + rates 




M,, 2 (GeV) M„ 2 (GeV) 



Figure 4: A\ 2 contours showing Mq against M1/2 for 14 TeV/1 fb _1 data. Fits are based on the 
four standard edges of group I with rates (upper right), and on the observables of groups I and II 
with (lower right) and without rates (lower left). 

the stransverse mass, does not lead to a significant improvement. The rate information 
together with the standard edges thus seems sufficient to determine the parameters of an 
mSUGRA fit to SPSla-type SUSY scenarios in the initial phase of a high energy LHC run. 
Our results for 14 TeV and 1 fb _1 without rates, presented in Table 2 and Figure 4, are 
consistent with the results presented in Ref. [15], specifically Table 15 thereof. 

At 14 TeV and 10 fb _1 also the systematical errors from lepton and jet energy scales 
are expected to improve, leading to a further reduction of the error on edges. Fitting the 
four edges of group I and the rates leads to a very accurate determination of the mSUGRA 
parameters Mq and -Mi/ 2 and of tan f3. At high energy and high luminosity, the observables 
of group III may become accessible. With all 13 observables of groups I, II and III included, 
the fit to the SPSla mSUGRA scenario is so well constrained that adding rates does not 
lead to a significant improvement, see Figure 5. Our results are consistent with those 
presented in Table 16 of [15]. 

3.3 Non-universal mSUGRA fits 

Let us now briefly discuss a more general class of fits where we attempt to determine the 
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SPSla 


M [GeV] 
100 


Mi [GeV] 
250 


M 2 [GeV] 
250 


M 3 [GeV] tan/3 
250 10 


A [GeV] 
-100 


7 TeV and 1 fb" 1 

I + rates 


Q1 1 + 2 7-3 
yi - i -36.1 


OQ« K +67.1 

236.5 I 579 


242 

— OO. I 


251.0±1;| 10.5 + _H 


-6.0 tl°sf 



14 TeV and 1 fb _1 

I + rates 
I + II, .rates 
I + II + rates 



98.5 till 
102.7 

98.6 ±i?;| 



245.8 ^ 
258.0 

249.6 + _f£ 



244.2 

255.4 ±£* 
248.7 



250.3 
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-7.0 



251.4 1? 2 9 2 



+6.0 
-7.1 



6.2 
10.3 



+11.0 
-2.0 
+5.9 
-3.0 



252.1 



9.5 +_r 7 



-389.9 tllf 
-102.0 tZ 7 

-12 7 5 + 790 
izi.o _ 204 



14 TeV and 10 fb -1 

I + rates 
I + II, xates 
I + II + rates 
I + II + III, urates 
I + II + III + rates 



qq 9 +13.3 
—9.3 



ocq 7 +37.1 

zoa. c _ 3a5 



+7.9 
-4.1 



98.4.0 
102.0 t 2 S 
99.1 11'° 

99.1 ±g : f 



246.8 
258.3 
245.9 
251.3 



+22.5 

-13.0 

+9.2 

-20.0 

+19.5 

-8.9 

+17.5 

-11.2 



256.9 
248.0 
254.3 
248.4 
250.4 



+20.2 
-23.9 
+12.6 
-6.8 

+6.0 
-10.2 

+10.3 
-5.2 

+9.9 
-5.1 



250.9 

249.2 
251.9 
248.7 
251.2 



+8.6 
-5.7 
+5.2 
-3.2 
+2.1 
-5.4 
+5.2 
-2.2 

+3.0 
-4.6 



14.3 
9.6 
9.9 
9.9 
9.9 



+1.7 
-9.7 
+1.7 
-0.8 
+1.4 
-1.0 
+0.9 
-0.7 

+1.1 
-0.7 



186.6 t^ll 



-117.5 



+83.1 
-45.5 



-124 7 + 101 

-63.6 



-98.6 
-101.0 



+41.7 
-48.0 
+47.9 
-50.5 



Table 3: Fits to mSUGRA parameters with non-universal gaugino masses for SPSla, with 
("+ratcs") and without ("xateS") event rates as an observable. I, II and III refer to the inclusion 
of the groups of previously considered observables (mainly edges) defined in the main text. 

GUT-scale gaugino mass parameters Mi , M 2 and M3 as fit parameters individually, instead 
of a common mass M 1 / 2 - Models without gaugino mass unification at the GUT scale can 
lead to very different phenomenology, see e.g. [62, 63], and have, to our knowledge, so far 
not been considered in SUSY parameter fits of LHC data alone. As in section 3.2, we 
use the SPSla observables collected in Table 1 as input. Even though SPSla is defined by 
Mi = M 2 = M3 = 250 GeV at the GUT scale, our results presented below demonstrate the 
importance of rates for the determination of individual gaugino mass parameters, crucial 
for the analysis of more general, non-universal, models. 

Analogous to what we observe for the fit of universal mSUGRA parameters, we find 
that a measurement of the four standard edges of group I alone is not sufficient to reliably 
determine the SUSY parameters of non-universal models in the initial phase of LHC at 
7 TeV and with 1 fb _1 . This is born out by the profile likelihood contours shown in 
Figure 6 (left). Adding rates, however, we obtain a reasonable fit of the larger parameter 
space including M±, M 2 and M3, as shown in Table 3 and Figure 6 (right). We find an 
accuracy of about 25% on Mi, M 2 and 5% on M3, see Table 3. The high accuracy on M3 is 
a result of the sensitivity of the production cross sections to the gluino mass, cf. Figure 2. 

Going from 1 fb" 1 at 7 TeV to 1 fb" 1 at 14 TeV leads to a significant reduction of 
uncertainties. With rates and the observables of groups I and II included, we obtain a good 
fit of the non-universal model, with errors in the range of 10% for Mq, Mi and M 2 , and 
about 3% for M3. The fit to tan/3 has an accuracy of 50%, compared to 30% for the case 
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of universal mSUGRA fits. With 10 fb" 1 at 14 TeV, finally, one can determine M3 with 
an error of less then 2%, while the uncertainties on Mq, Mi and Mi are in the range of 
5%. Considering the observables of groups I and II together with event rates is sufficient 
to arrive at the final LHC accuracy for the GUT-scale gaugino mass parameters Mi, M2 
and M3, see Table 3 and Figure 8. 
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Figure 5: A\ 2 contours showing Mq against M\/2 for 14 TeV/10 fb^ 1 data. Fits are based on 
the four standard edges of group I with rates (upper right), on the observables of groups I and II 
with (middle right) and without rates (middle left), and on observables of groups I, II, and III with 
(lower right) and without rates (lower left). 
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Figure 7: A\ 2 contours showing Mq against M3 for 14 TeV/1 fb" 1 data. Fits are based on the 
four standard edges of group I with rates (upper right), and on the observables of groups I and II 
with (lower right) and without rates (lower left). 
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Figure 8: A^ 2 contours showing Mq against M3 for 14 TeV/10 fb _1 data. Fits are based on 
the four standard edges of group I with rates (upper right), on the observables of groups I and II 
with (middle right) and without rates (middle left), and on observables of groups I, II, and III with 
(lower right) and without rates (lower left). 
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4. Conclusions and outlook 

Cross sections and branching ratios provide important information on the Lagrangian pa- 
rameters of TeV-scale super symmetry. We have presented a new method to include event 
rates, i.e. cross sections x branching ratios x cut acceptances, into global fits of SUSY pa- 
rameters at the LHC. While we expect event rates to be particularly important for SUSY 
scenarios that are not well constrained by the measurement of the standard kinematic 
edges in cascade decays, we have demonstrated that cross sections and branching ratios 
also add important information to fits of SPSla-type minimal supergravity models. In par- 
ticular, we find that event rates are crucial for a reliable determination of the mSUGRA 
parameters in the initial phase of LHC data taking at 7 TeV collision energy with 1 fb _1 
integrated luminosity. We have also studied a more general class of models where we allow 
for individual, in general non-universal, gaugino masses Mi, M2 and M3 as fit parameters 
instead of a common mass M]y 2 - The parameter determination of such models improves 
significantly when rates are included. 

The purpose of this paper is to establish a new method for including event rates into 
global SUSY parameter fits at the LHC and to quantitatively study their impact using 
SPSla as a case study. Our numerical results are based on Fittino, and the calculation 
of event rates will be included in the next release of the Fittino program package. Note, 
however, that the method is general and can also be used with other fitting codes. There 
are many possible studies that could follow from the method and the results presented here. 
We shall, for example, extend the range of applicability of our cut acceptance calculation by 
including three-body decay modes and decays into on-shell Z and Higgs bosons. Moreover, 
the estimated uncertainty on the acceptance should be verified by a detector simulation. So 
far, we have focused on LHC data alone, but shall address the interplay of LHC results with 
low-energy observables, other collider data and dark matter constraints in future analyses 
(cf. [13, 15, 64, 65, 66, 67, 68]). We have seen that event rates are particularly important 
to analyze more general models, such as models with non-universal gaugino masses. We 
thus plan a more comprehensive study of the determination of the gaugino mass pattern 
from LHC data, which may play a crucial role to distinguish models of supersymmetry 
breaking. 
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